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Abstract 

First time-resolved photoemission experiments employing attosecond streaking of electrons emit- 
ted by an XUV pump pulse and probed by a few-cycle NIR pulse found a time delay of about 100 
attoseconds between photoelectrons from the conduction band and those from the 4f core level of 
tungsten. We present a microscopic simulation of the emission time and energy spectra employing 
a classical transport theory. Emission spectra and streaking images are well reproduced. Differ- 
ent contributions to the delayed emission of core electrons are identified: larger emission depth, 
slowing down by inelastic scattering processes, and possibly, energy dependent deviations from the 
free-electron dispersion. We find delay times near the lower bound of the experimental data. 

PACS numbers: 79.60.Bm,72.10.-d,78.47.J- 
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I. INTRODUCTION 



Photoemission from solid surfaces employing photon energies ranging from the extended 
ultraviolet (XUV) to X-ray energies is a well established diagnostics tool to explore tjhe 
electronic band structure [I], the crystal structure 2|], and its chemical composition 



a. 



High-resolution angle resolved photoemission can provide detailed information on band gap 



S allows probes of the local 
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5|]. Until recently, virtually 



and conical intersections, e.g. in single-layer graphene [4]. EXA 
environment by spectral fluctuations due to interfering paths 
all of these techniques were based on the extraction of spectral information only, i.e. on the 
energy distribution -P(-E'e) for a given initial photon energy E^. The use of pulsed sources on 
the femtosecond scale has recently given access to dynamical processes, e.g. the formation 
of image states [g]. ^ 

The introduction of the attosecond streaking technique originally developed for gas- 
phase photoionization, to photoemission from solid surfaces has opened up a new perspective: 
observation of electronic motion in condensed matter and near surfaces in real time. The 
first proof of principle experiment revealed that XUV photoemission from a tungsten surface 
features an intriguing time structure with the conduction band electrons coming first while 
electrons from core levels, in the present case mostly 4f and 5p, are delayed by about 100 
attoseconds (110 ± 70 as) fl]. 

Analyzing the energy-time spectrum and identifying possile sources of the delay poses a 
considerable challenge to a microscopic simulation: apart from the intrinsic difficulties in 
accounting for the many-body response on a sub-femtosecond scale, the near infrared (NIR) 
probe with intensities of / = 2 ■ 10^^ W/cm^ is sufficiently strong as to actively modify the 
electronic response and emission process. Deviations from a field-free emission scenario are 
therefore to be expected. We present in the following a simulation of the attosecond streaking 

n 

of tungsten employing a classical transport theory (CTT, [9']). It allows to account for elastic 
and inelastic scattering processes as well as modifications of emission energy and angular 
distributions due to the presence of the streaking field. We find, overall, good agreement 
with the attosecond-streaking emision spectrum. The time spectrum shows, indeed, delayed 
photoemission with a delay time r ^ 40 as being at the lower bound of the experimental 
value. We present a detailed breakdown of different processes taken into account influencing 
the time spectrum and discuss possible additional contributions. 
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II. SIMULATION 



Before presenting the key ingredients entering the theoretical treatment, we briefly review 
the experiment underlying the scenario of the present simulation j^. 

Two coUinear linearly polarized laser pulses with the polarization direction in the plane 
of incidence were directed on a W(llO) surface under a grazing angle of incidence (Fig. 
[1]). The two pulses were an XUV pulse at 91 eV with a full width at half maximum 
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FIG. 1: (Color online) Schematic view of experimental setup: A 300 as XUV pulse hits the surface 
under a gracing angle of incidence Bin — 75.5° ionizing target atoms. Excited electrons propagate 
in the field of a 10 fs NIR pulse and are detected by a TOF spectrometer mounted perpendicular 
to the surface with an acceptance cone of half-width A9 = 5°. 



(FWHM) of ~ 6 eV and a pulse duration txuv — 300 as and a NIR pulse at 702 nm 
with ttv/j? — 10 fs. The angle of incidence was chosen to be at the Brewster angle for the 
probe pulse in order to avoid the influence of the deflected pulse on electrons emitted from 

n 

the surface (6'm = arctan(nAr//j) ~ 75.5°; n^jR = 3.85, kNiR = 2.86, [10|; dielectric function 



e = 81+182 



[n 



'k'^) + 2ink). This implies for the XUV pulse {nxuv = 0.93, kxuv = 0.04; 



10l |) total reflection at the surface and excitation of photoelectrons by an evanescent wave. 
Electrons excited by the pump pulse propagate in the field of the probe pulse modulating 
their energy as a function of their emission time ("streaking", ^]). Photoelectrons escaping 
the target surface were detected by a time-of-flight spectrometer with the detection direction 
normal to the surface. 40 streaking spectra were recorded with relative delays between pump 
and probe pulses ranging from —6 fs to +5 fs with respect to the maxima of the envelope 
functions of XUV and NIR pulses. 

On top of a strong background signal originating from above-threshold ionization (ATI) 
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two prominent features were observed in each of the streaking spectra: a lower energy peak 
around 55 eV and a higher energy peak at 85 eV. The former was attributed to electrons 
excited from 4f core states of tungsten, the latter to electrons from the target conduction 
band (consisting of 5d and 6s states). 

The incident XUV spectrum has been determined by measuring the kinetic energy of 
electrons emitted from a Ne gas target (Fig. [2]). Two main peaks from emission of 2s and 2p 
electrons with binding energies of 48.5 eV and 21.65 eV, respectively, can be distinguished. 
We have fitted the peak at 70 eV by a Gaussian function (solid line) from which the XUV 
photon energy of about 91 eV and the approximate FWHM of 6 eV were derived. 
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FIG. 2: (Color online) Kinetic energy of electrons emitted from Ne irradiated by the XUV pulse 
used in the present surface experiment (red dotted line); spectral distribution used in simulation 
(blue solid line). From the binding energy of the Ne(2p) level the XUV photon energy of ~ 91 eV 
was determined. 



The simulation involves the solution of the Langevin equation of motion, 

P= -FNmir, t) + Fstoc{t) + Fsurfir) 



(1) 



for an ensemble of a large number (typically 10^) of initial conditions representing the 
primary photoionization electrons. Fstoc{i) allows for both elastic and inelastic scattering of 
the liberated electron with tungsten cores as well as conduction electrons 



(2) 
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where Apj is the random colhsional momentum transfer as deduced from differential scatter- 
ing cross sections, ti are the random colhsion times determined from a Poisson distribution 
with a mean value given by the mean free path (MFP) or, equivalently, the mean flight 
time. The key point of this method is that the stochastic sequence (Apj,tj), determined 
either from quantum calculations or independent experimental data, allows quantum scat- 
tering information including diffraction to be included in an otherwise classical calculation 
(Eq. [H j^). In between collisions, free electrons with an effective mass (the value of which 
is discussed below) propagate within the penetration depth of the time-dependent NIR elec- 
tric field. Near the surface, the electrons must overcome the effective surface potential. 



the asymptotic part of which is the dynamical image potential This gives rise to the 
additional force -Fsur/(^) in Eq- [H 

As collisions can result in secondary electron emission, the simulation does not preserve 
the number of particles but allows "en route" for the generation of additional trajecto- 
ries thereby simulating the collision cascade. However, slow secondary electrons generated 
outside the energy window of interest are omitted in the following. 

A. Primary photoemission as initial condition 

The ensemble of initial conditions represents photoemission of tungsten electrons by a 91 
eV XUV pulse which originate from 5 atomic levels: 6s, 5d, 4f, 5p, and 5s. Photoelectrons 
from the latter (i?fej„^(5s) = 75.6 eV) are submerged in the ATI background and are therefore 
omitted in our simulation. Electrons excited from the 4f and 5p levels (referred to in the 
following as core levels) give rise to a peak near 55 eV kinetic energy. The binding energies 
of the 4f5/2 and 4f7/2 states are 33.5 eV and 31.4 eV, respectively. The branching ratio 
for photon energies around 100 eV of Rj = h/2/h/2 = 1-56 was taken from experimental 



photoemission data [l2l]. With a binding energy of 36.8 eV the energy of excited 5p3/2 
electrons overlaps with the distributions of 4f electrons. Due to the smaller photoexcitation 
cross section its relative importance is about a factor 3 smaller than the combined 4f levels 
[13I . Excitation of the 4f and 5p levels with the XUV pump pulse taking the 4f-branching 
ratio into account leads to a slightly asymmetric energy distribution of primary electrons 
around 58 eV. For brevity, we refer to this overlapping distribution in the following as the 
core level. 
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The 5d (4 electrons per atom) and 6s (2 electrons per atom) levels of tungsten form the 
conduction band. Its density of states (DOS) was taken from ij] (Fig. [3l) . The Fermi energy 
is 9.75 eV while the workfunction of W(llO) is taken to be 5.25 eV 15|. We split the DOS in 
two components: a free-electron gas like subband containing the 6s electrons (DOS oc \fE\ 
area below red line) and a more localized component representing the 5d electrons. A similar 
decomposition of the DOS with a slightly higher value for the Fermi energy has already been 



proposed earlier by Mattheiss 



16| . He suggested the bottom of the 6s band to lie below the 



bottom of the 5d band, as shown in Fig. [31 Density-functional theory calculations using the 
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FIG. 3: (Color online) Occupied part of the density of states of tungsten [l^. States below the red 
solid line (shaded area) originate from the 6s band, states above the line from the (more localized) 
5d level. 



18j give a Fermi energy 



ABINIT package [17| with pseudopotentials for the tungsten cores 
of Ep = 9.6 eV in good agreement with Q and Q. 

We note that in the spectral region of interest an additional emission channel may con- 
tribute: Auger decay from the conduction band filling a 5s hole created by photoionization 
(binding energy of the 5s level ~ 75.6 eV). The expected energy range for this Auger emis- 
sion would be 45 to 65 eV. We expect, however, that the refilling of 5s holes will proceed 
predominantly via the 5s5p^ decay channel. This is also supported by XPS spectra of the 
W 4f region in which no signatures of Auger electron emission have been seen (e.g. jigl). 
Therefore, we do not include this channel in our simulation. 

We set the emission strength (dipole oscillator strength) for each level to be constant 
within the spectral width of the XUV pulse. For the relative intensity ratios between differ- 
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ent levels, we have explored two different options. One uses the calculated photoionization 
cross sections for atomic tungsten [13]. Accordingly, the photoionization cross section for 
the 4f subshell of atomic tungsten is by about a factor two larger than for the 5d level and 
by more than an order of magnitude larger than for 6s electrons. Using this input we find a 
spectral distribution (in absence of any streaking field) in strong disagreement with experi- 
mental data. Experimental photoemission cross sections from the conduction band appear 
to be considerably larger than estimated from atomic cross sections. Alternatively, we have 
adjusted the strength of the conduction band peak relative to the core peak in order to 
reproduce a photoelectron spectrum after excitation by an XUV pulse (in the absence of the 
NIR field - thick red line in Fig. Hj). As the experimental spectrum includes the modification 




energy [eV] 

FIG. 4: (Color online) Energy spectra of electrons emitted by an XUV laser pulse in absence of 
the NIR streaking pulse from a tungsten surface. Experimental data with the ATI background 
subtracted are indicated by thick solid lines (blue - raw spectrum, red - smoothed raw spectrum), 
results of the CTT simulation by the thick dashed line. The thin dashed line represents the 
excitation spectrum as represented by the initial conditions. 

of the primary spectrum due to multiple scattering during transport until electrons cross the 
metal-vacuum boundary, a proper comparison must account for these processes as well (for 
details see below). Including transport and broadening of the emission spectrum to simulate 
the detector resolution (Gaussian distribution, a = 2.5 eV) the simulated photoemission 
spectrum (thick dashed line) fits the experimental data (raw data with the ATI background 
subtracted) rather well. A notable exception is the low-energy shoulder of the 80 eV con- 
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duction band peak the origin of which is presently not known. We note parenthetically that 



20|. 



a very recent measurement shows a much reduced shoulder 

The angular distribution of primary photoelectrons excited by linearly polarized photon 
beams is given in dipole approximation by 

= ^ (1 + /3(^)P.(co= 9)) = (l + ^(3 co.^ 9-1)) (3) 

where (T„^ is the total photoionization cross section from the ni subshell, P2 is the sec- 
ond order Legendre polynomial, and i3{uj) is the energy- dependent asymmetry parameter. 
(3 = describes isotropic emission while a value of /5 = 2 gives a pure cosine distribution. 
Calculated values for beta for the considered levels vary from (3 = 0.7 (5p) to /3 = 2 for 



s- levels 



13| . Other sources, however, suggest different values and experimental data are, to 



our knowledge, missing. We have therefore performed simulations for the limiting cases of 
(3 = and (3 = 2 but have found no significant influence of the specific choice on our final 
results. As an evanescent wave propagates parallel to the surface, its polarization vector is 
perpendicular to the surface (9 = 0). All results presented in Sec. [TTllwere calculated using 
(3 = 2 (cosine distributions along surface normal) for all levels. 

For the modification of the primary photoionization spectrum by multiple scattering and, 
in turn, for analyzing the time structure underlying the emission spectrum, the spatial distri- 
bution of the primary excitation, i.e. the depth profile of the source is of crucial importance. 
Starting positions of electron trajectories were distributed over the 20 topmost layers of a 
W(llO) crystal (lattice parameter aw = 3.16 A, layer spacing in (110) direction d = 2.24 A). 
For the free-electron like 6s component of the conduction band (Fig. [3]) we use the jellium 
approximation, i.e. constant density within the target material starting half a layer spacing 
above the topmost atomic layer. For the 5d subband of the conduction band as well as 
for the core electrons localized source distributions were chosen from Gaussian distributions 
with a FWHM of half a layer distance centered at crystal layer positions in agreement with 
our ab-initio calculations. 

The probability of excitation was assumed to be proprtional to the attenuated intensity of 
the XUV laser light entering the crystal. For 6in = 75.44° and nxuv = 0.93 total reflection 
conditions are fulfilled and only an evanescent wave enters the target. The decay length 
is about 6 ~ 0.55 ■ A ~ 75.6 A or, equivalently, more than 30 layer spacings. The depth 
dependent reduction in ionization probability has been taken into account, the first five 
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layers of which are shown in Fig. [51 
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FIG. 5: (Color online) Depth-dependent spectral distribution of primary photoelectrons excited 
by the XUV pulse with mean energy 91 eV, shown for the five topmost crystal layers. Projections 
onto the energy axis and the depth axis indicate the spectral and depth profiles. 



Electron transport 

Electron transport within the target material is calculated using a classical transport 



simulation described in detail in |2ll]. In brief, an electron released at its starting point 
with initial kinetic energy Ekin is subject to elastic and inelastic scattering processes as well 
as to deflection of its trajectory due to the influence of the electric field of the probe laser 
pulse in the target (Eq. [1]). The latter is refracted and damped by the target material. For 
fiNiFt. = 3.85 and k^jji = 2.86 at A = 702 nm we find a propagation angle of 6' ^ 14.5° with 
respect to the surface normal. The NIR intensity after a propagation length I is described 
by the Beer-Lambert-Law 

/(/) = /o ■ exp . /| = Jo . exp(-a/) (4) 

with the damping constant a = 0.005 corresponding to a penetration depth of about 
85 layers into the target. 

Elastic scattering cross sections have been calculated with the ELSEPA package 221] using 



a muffin-tin potential for the crystal atoms. From the energy-dependent total cross section 
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and the density of tungsten atoms the elastic mean free path (EMFP) is derived (Fig. [5]). 
The scattering angle in an elastic scattering event is determined by the energy- dependent 
differential cross sections. 

The inelastic scattering mean free path (IMFP) is derived from the momentum and energy 
dependent dielectric constant of the material e(g, ui) as j^ ] 



1 



dqduo TiEq 



1 



e{q,uj) 



(5) 



while the angular distribution of inelastically scattered electrons follows from 24 1 

1 



du! 

12 



UJ 



Q[E-E 



f 



(6) 



d^l _/ g2 Y ^ \ £{q,uj 

In Eqs. |5]and[6]i? is the instantaneous energy of the electron measured relative to the bottom 
of the conduction band while u and q are the collisional energy and momentum transfers, 
respectively. The B (step) functions impose the constraints of energy and momentum con- 
servation in the scattering event. e{q,uj) is constructed from an extrapolation of the optical 
data [e{q = 0, u)] for tungsten |loj] to the q — u plane (e.g., 25, 26]). The resulting inelastic 
mean free path (Fig. [6]) in the energy range from 50 to 100 eV is larger than the EMFP by 
factors 3 to 10. 
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FIG. 6: Calculated elastic and inelastic mean free paths of electrons in tungsten. The energy is 
measured from the bottom of the conduction band {Ep + IF ~ 15 eV). 



Due to the smaller EMFP the average escape depth for electrons emitted without having 
suffered energy losses will be smaller than the IMFP. This enhances the surface sensitivity 
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of the photoemission process. The total mean free path is determined by 

III 

■^tot ^el '^inel 

with \tot/ \nei,ei being the probabihty for an (in-)elastic scattering process to happen. Be- 
tween subsequent scattering processes electrons propagate in the time-varying electric field 
of the NIR laser. Travelling time, energy, and direction of motion are constantly updated. 
If an electron reaches the surface of the target it has to overcome the surface barrier which 
leads to an additional deflection at the surface. Energy lost in an inelastic scattering event 
is transferred to a secondary electron originating at the position of the collision. Its starting 
time is given by the time elapsed between start of the primary electron and the scattering 
event. As soon as the energy of an electron drops to below 25 eV the trajectory calculation 
is stopped. 

The effective mass of the energetic electron propagating through the crystal or, equiv- 
alently, its group velocity Vg as determined from the dispersion relation E{k), may have 
an important influence on the time dependence of the photoemission. To investigate the 
possible influence of the dispersion relation on the observed delay of emitted electrons we 
consider two limiting cases: a) a free-particle dispersion relation, E = k'^/2, with effective 
mass rrieff = 1 a.u. (red line in Fig. [T]) and b) the distribution of group velocities (green line 
in Fig. [7]) along the (110) direction of W calculated by Silkin et al. (supplementary material 
to ^). The latter features dips in the group velocity distribution around 65 eV and 120 eV 
with respect to the bottom of the conduction band due to the crystal potential. At both 
the expected energies of the core electrons and the conduction band (when measured from 
the bottom of the conduction band) we find Vg very close to that of a free electron (Fig. [7j). 
In Ref. js], however, it was assumed that the energy of the core electrons would coincide 
with the local minimum of Vg near 65 eV (see Fig. [Tj). Under this assumption larger run 
time differences between core electrons and conduction band electrons would result. For test 
purposes and in order to maximize the effect of the crystal dispersion relation we followed 
this prescription and shifted the envelope of group velocities to match its local minimum 
at 65 eV with the energy of the core electrons near 73 eV. The largest run-time differences 
reported in Sec. IIIII were found under these somewhat arbitrary assumptions. 

Finally, electrons escaping the surface are subject to the streaking field of the NIR laser 
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FIG. 7: (Color online) Group velocity distributions used in simulation; Vg = k {m^ff = 1) (red 
line), envelope of calculated group velocities (green line) from the dispersion relation of W along 
the (110) direction. Energies are measured with respect to the bottom of the conduction band; 
Ep + ~ 15 eV with Fermienergy Ep and workfunction W . 

pulse. The latter transfers a momentum of 

Ap= E{t)dt (8) 

to the electron where the integral is taken from the time of excitation to the conclusion 
of the laser pulse. Electrons with a final momentum oriented perpendicular to the surface 
with an acceptance angle AO = ±5° about the surface normal are included in the streaking 
image. Parameters describing the NIR pulse (maximum field strength, wavelength, pulse 
duration) were taken from 



III. RESULTS 



A. Emission energy and emission time 

We first present results for energy and emission time distributions in the absence of 
the streaking field. Energy and emission time spectra are strongly affected by energy loss 
properties during transport to the surface. Fig. M depicts the two-dimensional correlated 
primary excitation - final escape energy distribution. In the absence of inelastic processes 
the electron distribution would be located on the diagonal. Due to collisions described by 
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FIG. 8: (Color online) Primary excitation energy vs. escape energy for photoemission by an XUV 
pulse {huj = 91 eV) from a W surface. The loss function of tungsten (right-hand side) shows a 
pronounced excitation (loss channel) at about 25 eV. 

the energy loss function with a pronounced peak at about 25 eV, a significant portion of 
the initial distribution escapes at lower energy. Accidentally, the energy loss peak closely 
matches the energy spacing between core and conduction band electrons. Therefore, the 
primary conduction band electrons having undergone a single inelastic collision overlap with 
the primary energy distribution of the core electrons. The average excitation depth and 
traveling time in the target material for such electrons can be larger and will therefore 
contribute to the observed time delay. Also the secondary-electron background in the escape 
spectrum will be much larger in the spectral region of the core electrons adding to the 
observed run-time difference. 

The two-dimensional correlated escape-time escape-energy distribution (Fig. [9]) after the 
excitation by the XUV pulse feature consequently a much broader and slightly shifted escape 
time distribution of the "core" electron peak near 58 eV compared to the conduction-band 
distribution near 83 eV. Note that the tesc = line corresponds to the temporal maximum 
of the XUV amplitude. The maxima of the distributions at 58 and 83 eV are located close 
to 115 and 100 as, respectively. Averaging over the entire peak areas (44-70 eV, 70-95 eV) 
results in run times of 157 as and 115 as (see horizontal lines in Fig. [TOl) . The run-time 
difference (~ 42 as) can be accounted for in part by a simple estimate using the material 
properties of W: the average escape depth of electrons is given by their IMFP (see Fig. |6]). 
Using the free-electron dispersion relation the total travel time of electrons with energies 
of 73 and 98 eV (measured from bottom of conduction band) along their respective IMFP 
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FIG. 9: (Color online) Correlated escape-time escape-energy distribution after excitation of a W 
surface by an XUV pulse {hw = 91 eV). tesc = corresponds to the peak amplitude of the XUV 
pulse. 
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FIG. 10: (Color online) Energy spectrum (thin line) and escape times of electrons (thick red line) 
excited by an ultrashort XUV pulse. Horizontal lines indicate the escape time averaged over the 
core and conduction band peak regions. The mean run-time difference is about 40 as. 

of 6.5 A and 5.8 A is about 125 as and 100 as, respectively. This gives a lower bound of 
the run-time difference for electrons with escape energies of 58 eV and 83 eV of about 25 
as. This runtime difference is increased by scattering events in the target material which 
primarily increase the average run time of electrons ending in the spectral region of the core 
peak. 

Projection onto the energy axis provides the spectrum while temporal information (Fig. 

14 



[TOi) is extracted by averaging the escape time over the electrons within a given escape energy 
bin. While at the maximum of the energy spektrum the run-time difference is only about 20 
as, averaging over the entire energy range of the core and conduction band peaks (horizontal 
lines) results in a difference of 42 as. This is at the lower bound of the experimental value of 
110 ± 70 as. The agreement would be considerably better when using the shifted envelope 
of group velocities as discussed above. In this case the run-time difference would be more 
than double and would increase to 101 as. 

B. Streaking images and center-of-mass motion 

The time structure of the photoelectron emission is extracted in the experiment by 
attosecond streaking, i.e. by time-to-energy mapping in a few cycle NIR pulse. As soon as 
photoelectrons are excited they are subject to the NIR streaking field. If a free electron 
interacts with the entire laser pulse its momentum remains unchanged. If, however, it is 
set free while an electric field amplitude E{t) is present at the position of excitation, a 
net momentum is transfered (Eq. [8]). Varying the delay time between XUV pump and 
NIR probe pulses result in an oscillation of the final electron energy {t^ . In the experiment 
streaking images have been recorded in delay-time steps of At = 100 as. r < imply the 
probe pulse to precede the pump pulse, positive r signify the pump pulse preceding the 
probe. 

The characteristic streaking oscillations in the simulated electron energy can be clearly 
seen (Fig. [TTl) . From the streaking image a mean run-time difference of 33 as between core 
electrons (integrated over the energy interval < E < 70 eV) and the conduction electrons 
(70 < E <95 eV) can be deduced. The values deduced from the streaking image is slightly 
lower (by about 10 as) than that directly deduced from the simulated time spectrum. This 
discrepancy is not due to the resolution limits of the streaking technique but consequence 
of the modification of the electron transport by the streaking field. Electron trajectories are 
deflected in the direction of the NIR laser polarization inside the target. As a consequence, 
the escape depth of excited electrons is reduced. The point to be emphasized is that the 
streaking field not only probes but actively modifies the emission time spectrum. Along the 
same lines, the streaking field also alters the final emission energy spectrum (Fig. fT^ . Taking 
the streaking field into account, the shape of the electron spectra drastically changes. Almost 
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FIG. 11: (Color online) Calculated streaking images for electrons emitted by a 300 as XUV pulse 
and streaked by a few-cycle NIR pulse. Emission along the surface normal with an acceptance cone 
angle of 6* = 0±5°. 
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FIG. 12: (Color online) XUV pulse induced photoemission spectrum in the presence of a NIR 
streaking field. Thick red line: simulated spectrum for electrons emitted in the direction of the 
surface normal with an acceptance cone of half-width = 5°; black line: experiment; blue line: 
smoothed experimental data. 

perfect agreement between simulation and experiment is achieved. Both the simulation and 
the experiment feature broadened core and conduction band peaks (compare Figs. [5] and 
[T2I) . The overall agreement in peak height and width indicates that the present simulation 
is capable of accounting for most of the relevant processes governing the XUV-pulse induced 
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photoemission from a tungsten surface. 



IV. CONCLUSION 

We have presented simulations of the electron excitation, transport, and attosecond 
streaking for an XUV-pump-NIR-probe setting near a tungsten surface. This scenario mod- 
els the recent experiment by Cavalieri et al. |8| . We find excellent agreement for the energy 
spectrum in the presence of the streaking field. In agreement with the experiment we find 
the core electron emission to be delayed relative to the conduction band electrons. Calcu- 
lated run-time differences between groups of photoelectrons are at the lower bound of the 
error bar of the measurement. Latest experiments indicate a somewhat smaller run-time 
difference of 85 ± 45 as [20] reducing the gap between measured data and our simulations. 
Processes responsible for the time delay identified by the present simulation include the 
larger emission depth of core electrons, the contribution of primary emitted conduction elec- 
trons slowed down to energies matching emitted core electrons as well as secondary electrons. 
When we include a modification of the group velocity distribution for energies in the region 
of core electrons as proposed by Silkin et al. [8], the delay time increases and is closer to 
the experimental data. We note, however, that invoking this correction would require an 
energy shift of the spectrum for which a convincing explanation is missing. Possible other 
mechanisms not yet accounted for include the influence of the NIR pulse on the primary 
photoexcitation process and local crystal field effects on the emission time spectrum. The 
impact of both effects could be increased by local field enhancements at the surface (surface 
plasmons) as recently shown for the nanoplasmonic field microscope [27]. Such fields might 
influence excitation and transport of photoelectrons from the conduction band and bound 
states differently. 
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